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ABSTRACT 

We report on the discovery of gamma-ray pulsations from five millisecond pulsars 
(MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides pro- 
vided by various radio observatories. We also present confirmation of the gamma-ray 
pulsations from a sixth source, PSR J2051— 0827. Five of these six MSPs are in binary 
systems: PSRs J1713-I-0747, J1741-hl351, J1600-3053 and the two black widow bi- 
nary pulsars PSRs J0610-2100 and J2051-0827. The only isolated MSP is the nearby 
PSR J1024— 0719, which is also known to emit X-rays. We present X-ray observations 
in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly de- 
tected, we can only give upper limits for the X-ray fiux of the former. There are no 
dedicated X-ray observations available for the other 3 objects. 

The MSPs mentioned above, together with most of the MSPs detected by Fermi, 
are used to put together a sample of 30 gamma-ray MSPs which is used to study the 
morphology and phase connection of radio and gamma-ray pulse profiles. We show that 
MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission 
prcscmt the steepest radio spectra and the largest magnetic fields at the light cylinder 
among all MSPs. As well, we also observe a trend towards very low, or undetectable, 
radio linear polarisation levels. Those properties could be attributed to caustic radio 
emission produced at a range of different altitudes in the magnetosphere. We note 
that most of these characteristics are also observed in the Crab pulsar, the only other 
radio pulsar known to exhibit phase-aligned radio and gamma-ray emission. 

Key words: pulsars: general — gamma-rays: general — X-rays: general 



1 INTRODUCTION 

Pulsed gamma-ray emission from more than one hundred 
pulsars has been detected by the Large Area Telescope 
(LAT) aboard the Fermi Gamma-ray Space Telescope^ 
(Abdo et al. 2010b, 2013). About a third of them are mil- 
lisecond pulsars (MSPs), fast rotators that spin down very 
slowly and steadily. MSPs are thought to have been spun- 
up through accretion of material from an evolved compan- 
ion star (Alpar et al. 1982; Radhakrishnan & Srinivasan 
1982). During accretion, the system appears as a low-mass 
X-ray binary (LMXB) and no pulsed radio emission is ob- 
served. Although other possibilities have been considered 
(e.g. Ruderman 1991), it is commonly believed that accre- 
tion is responsible for the low surface magnetic fields ob- 
served amongst MSPs. In this model, the original field has 
boon buried by the in-falling material, decreasing its surface 
strength and decreasing the magnetic torque acting on the 
neutron star (Romani 1990; Bhattacharya & van den Heuvel 
1991). 

Most MSPs are found in binary systems. Among those 
in binaries, there is a growing number in very tight binary 
systems with low mass companions, the so called Black Wid- 
ows (e.g. Roberts 2011). It is believed that most solitary 
MSPs came through the black widow evolutionary path, in 
which the companion is ablated by the pulsar wind, cither 
during the X-ray phase or once the pulsar turns on in the 
radio (Eichler & Levinson 1988; Ruderman et al. 1989). 

The gamma-ray emission mechanism for MSPs and nor- 
mal pulsars is believed to be the same (Abdo et al. 2009). 
MSPs exhibit the lowest gamma-ray luminosities among 
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gamma-ray pulsars, a consequence of their low spin-down 
energy rates {E), compared to young pulsars (Abdo et al. 
2010b). Most observed gamma-ray pulse profiles consist of 
two dominant, sharp peaks, which suggest the emission is 
caustic in nature (Abdo et al. 2010b). Among the differ- 
ent models for radio and gamma-ray emission of pulsars, 
these properties tend to favour models in which gamma rays 
are generated in the outer magnetosphere (e.g. Romani & 
Yadigaroglu 1995; Dyks & Rudak 2003) and radio emission 
is produced at lower altitudes (e.g. Rankin 1993). This is 
supported by the fact that most gamma-ray detected ra- 
dio pulsars exhibit gamma-ray pulses arriving out of phase 
with the radio pulsations (Abdo et al. 2010b). However, this 
scheme has been challenged by the increasing number of 
MSPs found to exhibit phase-aligned radio and gamma-ray 
pulse profiles, suggesting that, at least in these cases, they 
are both produced at a similar location in the magneto- 
sphere (Abdo et al. 2010a; Freire et al. 2011a; Guillemot 
et al. 2012b). It has been noted that some of the MSPs 
showing phase-aligned radio and gamma-ray emission ex- 
hibit very low levels of radio linear polarisation (Guillemot 
et al. 2012b), which is predicted to some extent by caustic 
emission models (Venter et al. 2012). 

Before these discoveries, the Crab pulsar was the only 
pulsar known to exhibit phase-aligned radio and gamma-ray 
emission (Kuipcr et al. 2003). MSPs present radio polari- 
sation properties similar to normal radio pulsars (Xilouris 
et al. 1998; Yan et al. 2011; Keith et al. 2012) but often 
much wider pulse profiles. Interpulses and bridge emission 
between different peaks are more common among MSPs 
(Yan et al. 2011). MSPs present radio spectra with a slope 
similar to normal radio pulsars, although arguably slightly 
steeper (Toscano et al. 1998; Kramer et al. 1998; Maron 
et al. 2000). 

As well as detecting many new radio MSPs in gamma 
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rays (Abdo et al. 2009), Fermi has also pointed the way to 
many MSP discoveries via their gamma-ray emission prop- 
erties (e.g. Ransom et al. 2011; Cognard et al. 2011; Keith 
et al. 2011; Guillemot et al. 2012a; Kerr et al. 2012). Here we 
report on another six MSPs detected in gamma rays by the 
Fermi LAT, all of which were previously known radio pulsars 
(PSRs J0610-2100, J1024-0719, J1600-3053, J1713+0747, 
J1741-I-1351 and J2051-0827). This is the first time that 
pulsed gamma-rays are detected from these sources, except 
for PSR J2051— 0827, for which a 4cr detection was reported 
by Wu et al. (2012). 

By combining these six MSPs with 24 previously re- 
ported gamma-ray MSPs we are able to start evaluating 
trends in the emission properties of the gamma-ray MSP 
sample. In particular, we use this sample to study the pulse 
profile properties and other connections between the radio 
and gamma-ray characteristics of these objects. 

In section 2 we describe the methods used to detect 
and analyse the gamma-ray pulsations from the six new 
gamma-ray MSPs and in section 3 we comment on their 
multi-wavelength properties. In section 4 we present a search 
for orbital modulation in the gamma-ray emission of the five 
MSPs in binary systems. The study of the radio and gamma- 
ray properties of most known gamma-ray MSPs is presented 
in section 5 and the discussion of the study is in section 6. 
A summary of our work is in section 7. 



2 METHODOLOGY: DETECTING 
GAMMA-RAY PULSATIONS 

The detection of pulsed gamma-ray emission from known 
MSPs is possible through the use of precise rotational 
ephemerides obtained through frequent radio observations, 
which are used to assign a rotational phase to each gamma- 
ray photon (Smith et al. 2008). Photons are then binned in 
pulse phase to create a histogram which represents the light 
curve of the pulsar. 

2.1 Radio analysis 

We use pulse times of arrival (TOAs) obtained from radio 
observations with the Arecibo Observatory (AO, Preire et al. 
2011b), the Nangay Radio Telescope (NRT, Theureau et al. 
2005), the Jodrell Bank Observatory (JBO, Hobbs et al. 
2004), the Parkes Observatory (PKS, Weltevrede et al. 2010) 
and the Westerbork Synthesis Radio Telescope (WSRT, 
Karuppusamy et al. 2008; VoClte et al. 2002). See section 
2.1.1 for more details on the use of the different datasets. 

To study the phase alignment between gamma and ra- 
dio pulses it is necessary to account for the delay suffered 
by the radio waves in their passage through the interstellar 
medium. In order to do this precisely, accurate estimates 
of the dispersion measure (DM) are essential, especially for 
rapidly rotating pulsars (Smith et al. 2008) . With the aim of 
having accurate and updated DM values, we measured the 
DM using multi-frequency radio data taken during the time 
of the Fermi mission for five of the six gamma-ray MSPs 
(Table 1). 

The best available radio TOAs were used to produce 
radio ephemerides valid from at least 3 months before the 
start of the Fermi mission until October 2011. Ephemerides 



were produced by analysing the TOAs with the timing soft- 
ware Tempo2 (Hobbs et al. 2006). Because of the use of 
different pulse-profile templates to match the observations, 
hence different points of reference, and the existence of 
small clock divergences between observatories, the different 
datasets sometimes require phase-alignment. Thus, in ad- 
dition to the spin frequency and its first derivative, binary 
parameters, position, proper motion and DM, we also fitted 
for phase delays between the datasets from different obser- 
vatories. More information on the timing solutions built for 
each pulsar can be found in Table 1 and section 2.1.1. 

While pulsar timing can provide most of the known pa- 
rameters describing pulsars, distances can be obtained by 
different techniques. If no better distance estimate was avail- 
able (e.g. via parallax), the DM was used to calculate a dis- 
tance based on the Galactic Free Electron Model (NE2001, 
Cordes & Lazio 2002). The best distance estimates (d) for 
the 6 MSPs are listed in Table 2. We note that the error bars 
of the DM-based distances could potentially be largely un- 
derestimated, affecting the error bars quoted for gamma-ray 
luminosities and efficiencies (section 2.2.1). 

Because MSPs present low rotational period derivatives 
(P), the Shklovskii effect (Shklovskii 1970) can bo signifi- 
cant for those objects with a large proper motion (like PSR 
J1024— 0719, see below). We have corrected all P values for 
this effect and the corrected values can be found in Table 2. 
Other radio properties like position, rotational and orbital 
periods and proper motion are also found in Table 2. 

2.1.1 Building timing solutions, dispersion measures and 
distances 

For PSR J0610— 2100 most observations during the last 4 
years have been carried out at ~ 1.4 GHz, hence it is not 
possible to perform a sensitive analysis of the DM during the 
time of the Fermi mission. Nonetheless, using a few ~ 2 GHz 
TOAs and the small spread of observing frequencies around 
1.4 GHz of the other available TOAs, we are confident that 
the DM has not changed by more than 0.02% with respect to 
the value quoted in the radio discovery paper (Burgay et al. 
2006). This ensures an accuracy better than 0.006 rotations 
on the phases assigned to the gamma-ray photons. 

For PSR J1024-0719, we measure a DM of 6.488 ± 
0.001 pccm~^, which indicates a distance of 380 ± 40 pc, 
based on the NE2001 model. However, the Lutz-Kelker-bias- 
corrected distance is 500 ± 100 pc (Hotan et al. 2006; Ver- 
biest et al. 2012). As is discussed below, considerations of the 
transverse movement and the Shklovskii effect are important 
before adopting any of these values. There are slightly dif- 
ferent published proper motion values for this pulsar (see 
Bailes et al. 1997; Hotan et al. 2006). Using 5 years of NRT 
data at ~ 1.4 GHz we fitted for position and proper mo- 
tion and obtained values entirely consistent with those in 
Hotan et al. (2006) . The large total proper motion measured 
(~ 59 masyr^^) implies a significant contribution to the ob- 
served P value by the Shklovskii effect. In fact, it seems that 
most of the observed spin-down comes from this effect. The 
parallax distance indicates a negative intrinsic P and the 
maximum possible distance, in order to obtain a positive P, 
is 410 pc. Using the lowest possible distance given by the 
parallax uncertainties (i.e. 400 pc), we calculate the upper 
limit P ^ 0.05 x 10"^°. We wiU use this P upper limit and 
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Table 1. Properties of the radio ephemerides used to fold the gamma-ray data. 



Pulsar 


Observatories 


■'VtOAs 


RMS 


MJD 


range 


DM 


S4> 








fJ,S 






pccm-3 


10-3 


J0610-2100 


JBO, NRT, PKS 


111 


2.84 


54509 


- 55850 


60.67(l)t 


6 


J1024-0719 


JBO, NRT, WSRT 


145 


1.82 


54590 


- 55837 


6.488(1) 


0.4 


J1600-3053 


NRT 


193 


1.20 


54564 


- 55797 


52.3218(4) 


0.2 


J1713+0747 


JBO, NRT, WSRT 


257 


0.74 


54501 


- 55804 


15.9929(1) 


0.1 


J1741-I-1351 


AO 


8192 


0.77 


52840 


- 55889 


24.2014(1) 


0.04 


J2051-0827 


JBO, WSRT 


1590 


14.52 


54338 


- 55880 


20.73673(2) 


0.01 



Note. — Columns contain the observatories involved in the observations, the total number of 
TOAs used (A^xoAs): the root mean square variation (RMS) of the timing residuals and the 
time range of validity for the timing solution. The last two columns show the DM values and 
the uncertainty of the phases assigned to the gamma-ray photons caused by the uncertainties on 
the DM, 5(j>. DM uncertainties on the last quoted digits are in parentheses. These uncertainties 
correspond to the formal errors given by TEMP02 and may under-represent the real uncertainties. 
The photon phase errors (Sep) caused by the uncertainty on the DM were calculated using these 
values, t Value taken from Burgay et al. (2006). 



adopt the Lutz-Kelker-bias-corrected distance calculated by 
Verbicst ct al. (2012). 

For PSR J1600— 3053, the timing solution was obtained 
using data at two main frequencies (1.4 and 2.1 GHz), that 
allowed us to calculate an accurate and up-to-date DM. We 
also included the first time derivative for the DM in the fit 
(i.e. it was allowed to vary linearly with time), measuring a 
value consistent with that in You et al. (2007). 

For PSR J1713-F0747, a distance of 0.9 ± 0.1 kpc is ob- 
tained from the NE2001 model using our DM estimate (Ta- 
ble 2). However, using 12 years of radio timing, Splaver et al. 
(2005) measured this pulsar's yearly parallax and calculated 
a distance of 1.1 ± 0.1 kpc, which is the value we adopt here. 

For PSR J1741-fl351, the timing solution was obtained 
with observations carried out at the AO using the L-wide re- 
ceiver, recording data between 1.1 and 1.6 GHz, which offers 
enough bandwidth to monitor DM variations. Based on the 
NE2001 model and the DM obtained from these observa- 
tions the distance to this pulsar would be 0.9 ±0.1 kpc. The 
same datasot contains the clear effects of a yearly parallax, 
suggesting a similar distance of ~ 1.08 ± 0.05 kpc (this is a 
preliminary value, Freire et al. 2013), which is the value we 
use. 

For PSR J2051— 0827, the orbital parameters are chang- 
ing rapidly: Stappcrs et al. (1998) reported the measurement 
of the first time derivative of the orbital period, Doroshenko 
et al. (2001) detected up to 3 time derivatives of the orbital 
period and Lazaridis et al. (2011) detected variations of the 
orbital period and the projected semi-major axis. Almost 4 
years of WSRT data at 1.4 GHz was used to fit for position, 
proper motion and binary parameters, using the ELLl bi- 
nary model (Lange et al. 2001) implemented in Tempo2, 
and including one derivative of the orbital period. A sec- 
ond WSRT dataset taken at 0.35 GHz and JBO TOAs at 
1.5 GHz were then used together with the 1.4 GHz data to 
fit for the DM and its first derivative, keeping all orbital and 
astrometric parameters fixed. 

These timing solutions will be made available through 
the Fermi Science Support Center^. 



^ http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 



2.2 Gamma-ray analysis 

To study the gamma-ray emission from these six MSPs we 

selected Fermi LAT data taken between 2008 August 4 and 
2011 August 4 using the Fermi Science Tools (STs)'\ We re- 
stricted the dataset to events with energies between 0.1 and 
100 GeV, reconstructed directions within 15° of the pulsar 
locations, zenith angles smaller than 100°, and belonging to 
the "source" class of the P7_V6 instrument response func- 
tions (IRFs). We rejected the data collected when the LAT 
rocking angle exceeded 52°, when the instrument was not 
operating in the science observation (or configuration) mode 
or when the data quality flag was not set as good. 



2. 2. 1 Spectral properties 

The gamma-ray spectral properties of the pulsars were de- 
termined using a binned maximum likelihood method, as 
implemented in the pyLikelihood Python module of the 
Fermi STs. This method fits a model representing the point 
sources in the selected Region Of Interest (ROI) and the 
diffuse emission to the data, and finds the best-fit parame- 
ters optimising the likelihood function describing the data. 
Our models included all sources from the 2FGL catalog 
(Nolan et al. 2012) found within 20° of the pulsar posi- 
tions. Parameters of sources within 8° of the pulsars were 
left free in the fit, while parameters of sources more than 
8° away were fixed at the values listed in the 2FGL cata- 
log. The diffuse Galactic emission was modelled using the 
gaL2yearp7v6-vO map cube. The residual instrument back- 
ground and the diffuse extragalactic emission were modelled 
using the iso-p7v6source template*. The normalisation of 
the diffuse components were left free in the fits. The first 
step of the spectral analysis involved modelling the pulsar 
spectra with exponentially cutoff power-law (ECPL) shapes 



3 http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/ovcrview. html 
* These diffuse models are available for download from the Fermi 
Science Support Center, see http://fermi.gsfc.nasa.gov/ssc/ 
data/ access/lat/BackgroundModels . html 



New Fermi MSPs and the radio / gamma-ray connection 5 



Table 2. Main pulsar parameters for the 6 detected MSPs. 



Pulsar 


I 


b 


P 


P 


Pb 




d 


^corr 




Refs. 




(deg) 


(deg) 


(ms) 


10-20 


days 


mas yr" ^ 


kpc 


10-20 


xlO^'^ ergs~^ 




J0610-2100 


227.8 


-18.2 


3.86 


1.24 


0.3 


18.2(2) 


4(l)t 


0.1(3) 


1(2) 


1, 2 


J1024-0719* 


2.51.7 


40.5 


5.16 


1.85 




59.9(2) 


0.5(1) 


< 0.05 


< 0.1 


1, 3 


J1600-3053 


344.1 


16.5 


3.60 


0.95 


14.4 


7.2(3) 


2(1) 


0.84(4) 


7.1(4) 


4, 3 


J1713+0747 


28.7.5 


25.2 


4.57 


0.85 


67.8 


6.30(1) 


1.1(1) 


0.805(4) 


3.33(2) 


5 


J1741+1351 


37.89 


21.6 


3.75 


3.02 


16.3 


11.71(1) 


1.08(5) 


2.89(1) 


21.68(4) 


6 


J2051-0827 


39.19 


-30.4 


4.51 


1.27 


0.1 


7.3(4) 


1.0(2) 


1.21(1) 


5.21(5) 


7, 2 



Note. — The first 8 columns contain pulsar names, positions in galactic coordinates, periods, period derivatives, 
orbital periods, proper motions and distances. The 9*'' and 10*'' columns give the values for the period derivative 
and spindown energy rate, both corrected for the Shklovskii effect (Pcorr and Ecorr). Uncertainties in the last 
quoted digits are in parentheses. References in the last column are for proper motion and distance: (1) This work; 
(2) NE2001 model (Cordes & Lazio 2002); (3) Verbiest et al. (2012) (4) Verbiest et al. (2009); (5) Splaver et al. 
(2005); (6) Preire et al. (preliminary, 2013) (7) Laaaridis et al. (2011). ^See a discussion about this value at the end 
of section 2.2.1. *See a discussion about the distance to this MSP and the Shklovskii corrections in sections 2.1.1 
and 2.2.1. 



of the form: 
dN „ / S 



exp 



E 



(1) 



where A'o is a normalisation factor, F is the photon index 
and Ec is the cutoff energy of the spectrum. The spectral pa- 
rameters for PSRs J0610-2100 and J1713+0747 found with 
this model are listed in Table 3, along with the derived inte- 
grated photon fluxes F and energy fluxes G above 0.1 GeV. 
These pulsar spectra were also fltted with a simple power- 
law model to test the validity of the ECPL model by com- 
paring the goodness of the fit for the two models. For both 
pulsars the power-law model was rejected with more than 
3(7 significance. For the other four pulsars, fitting the spec- 
tra with all three parameters of the ECPL spectral shape 
left free was unsuccessful and yielded unsatisfactory results 
indicative of bad convergence issues. To estimate the cutoff 
energy and the integrated energy fiux G for these pulsars, 
a second spectral fit was performed using a ECPL model 
with two free parameters only, fixing F to a value of 1.3, 
which is the average seen for the 32 strongest MSPs in the 
second catalog of LAT pulsars (Abdo et al. 2013, in prepa- 
ration). The results of these fits are also listed in Table 3. 
In this table, the first quoted uncertainties are statistical, 
and the second are systematic. For PSRs J0610— 2100 and 
J 1713+0747, the latter uncertainties were estimated by run- 
ning the fitting procedure using bracketing IRFs where the 
effective area was perturbed by the estimated uncertainties 
±10% at 0.1 GeV, ±5% near 0.5 GeV, and ±10% at 10 GeV, 
using linear extrapolations (in log space) in between (Acker- 
mann et al. 2012a). For the other four MSPs, the systematic 
errors caused by fixing the photon index to a nominal value 
of 1.3 probably dominate the errors due to uncertainties in 
the IRFs. To estimate the systematic uncertainties resulting 
from this choice, we fitted the data using photon indices of 
0.7 and 2, these values representing the cxtrema observed for 
32 strong gamma-ray MSPs (Abdo et al. 2013). The best-fit 
parameters for F = 0.7 and F = 2 provide the limits on 
the actual values listed in Table 3. We note that the flts 
obtained with F = 2 were inconsistent with the data in all 
cases, which resulted in very large values of Ec for PSRs 
J1600-3053 and J1741+1351 and, therefore, large system- 
atic errors for these objects. Finally, the gamma^ray lumi- 



nosities L-y = A'KfciGdP and the efficiencies 77 = L-^jE for 
the conversion of spin-down luminosity into gamma-ray ra- 
diation were calculated assuming a beaming factor of /s; = 1 
(see Walters et al. 2009, for more details) which is common 
for outer magnetospheric emission models (e.g. Venter et al. 
2009). These quantities are given in Table 3. 

Two objects present anomalous efficiencies. We obtain 
r? > 0.8 for PSR J1024-0719, determined by the maximum 
E allowed by the P value adopted in section 2.1.1, which is 
very low. However, if this MSP was at 350 pc, E would be 
larger and thus the efficiency would be smaller and closer to 
values commonly observed (see the discussion on this MSP's 
distance in section 2.1.1). We note that this problem could 
also be alleviated if fn was smaller than 1. By modelling the 
gamma-ray emission of a few MSPs, Venter et al. (2009) cal- 
culated /fi for different geometric configurations (magnetic 
inclination and line of sight). Although for known gamma- 
ray MSPs they find /n values close to 1, under some specific 
geometric configurations their calculations show that /n also 
could be very small. 

The other case is the large efficiency obtained for PSR 
J0610-2100 (»? ~ 10, Table 3). An efliciency above 1 is un- 
physical and to have 7? < 1 requires ^2/0 < 1 kpc2. As noted 
above, it is possible that the flux correction factor is not 
exactly /n = 1, but smaller. We also consider the possibility 
that the large DM-based distance (Table 2) is due to mate- 
rial along the line of sight, not modelled in NE2001. Infrared 
images acquired by WISE (4.6, 12, 22 microns) and by IRAS 
(25, 60, 100 microns) show pronounced nebulosity around 
PSR J0610— 2100. A distance of 1 kpc in this direction cor- 
responds to 15pccm-^ in NE2001. Typical cloud sizes and 
over-densities can accommodate this DM discrepancy. For a 
distance of 1 kpc (taking /n = 1), the gamma-ray luminosity 
would be as low as 0.8x 10^^ ergs"^, the Shklovskii-corrected 
spin-down power would be 6 x lO"^"^ ergs"^ and the efficiency 
would be r/ = 0.1. A smaller distance also improves the 
transverse velocity estimate. The measured proper motion 
(Table 2) and a distance of 4 kpc imply a transverse velocity 
Vt = 345 km s"^, a rather large value, given that the mean 
value for MSPs has been estimated to be ~ 90kms-^, with 
a dispersion of 20kms-^ (Lyne et al. 1998). With a distance 
of Ikpc the transverse velocity of PSR J0610— 2100 would 
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be ~ 85kms^^, very similar to the mean value for MSPs. 
However, we note that if /n < 1 the distance (thus the 
transverse velocity) could be larger. For example, if we im- 
pose the maximum value of Vt allowed from the dispersion 
(i.e. Vt = 110kms"\ Lyne et al. 1998) then d = 1.3kpc 
and fn < 0.6. 

2.2.2 Search for pulsations 

It has been shown that pulsation searches can be made more 
sensitive by weighting each event by the probability that 
it originates from the considered gamma-ray point source 
(Kerr 2011; Guillemot et al. 2012b). These event probabili- 
ties depend on the spectra of the sources of interest and of 
the other sources in the region. The best-fit spectral mod- 
els and the Fermi ST gtsrcprob were thus used to com- 
pute the event probabilities. The event arrival times were fi- 
nally phase-folded with the Fermi plug-in distributed in the 
Tempo2 pulsar timing package (Hobbs et al. 2006; Ray et al. 
2011). The weighted H-test statistics (Kerr 2011) obtained 
by reducing the regions of interest to 5° around the pulsars 
are listed in Table 3. These test statistics values all corre- 
spond to pulsation significances larger than 5cr. We there- 
fore have detected pulsed gamma-ray emission from PSRs 
J0610-2100, J1024-0719, J1600-3053, ,11713+0747 and 
J1741+1351 for the first time, and confirmed the marginal 
detection of PSR J2051-0827 presented by Wu et al. (2012). 

Probability-weighted gamma-ray light curves for the six 
MSPs are shown in Figures 1 to 6, along with radio and 
X-ray profiles when available. The gamma-ray background 
levels in these figures were obtained by summing the prob- 
abilities that the events are not due to the pulsar, as de- 
scribed in Guillemot et al. (2012b). Statistical error bars 
were calculated as \/^~wf, where Wi is the event proba- 
bility and i runs over events in a given phase bin (Pletsch 
et al. 2012). In Table 3 we list the radio-to-garrirria-ray lags S 
and gamma-ray peak separations A for pulsars with multi- 
ple gamma-ray peaks, where the positions of the gamma-ray 
peaks were determined by fitting the integrated light curves 
above 0.1 GeV with Lorentzian functions, and the positions 
of the radio peaks were defined as the maxima of the radio 
light curves. 



3 MULTI- WAVELENGTH PROPERTIES OF 

THE SIX MSPS 

Below we present some basic information on these MSPs and 
detailed descriptions of their pulse profiles, in gamma rays 
and radio. If X-ray observations were available, descriptions 
of their main X-ray properties are also given. 

3.1 PSR J0610-2100 

PSR J0610— 2100 is a 3.8 ms pulsar in a 6.9-hour orbit with 
a ~ 0.02 M© white dwarf (Burgay et al. 2006). Given the 
low mass of the companion and the short orbital period, 
PSR J0610— 2100 is one of the Black Widow binary systems. 
However, in this case there are no radio eclipses nor DM 
variations caused by ablation of the companion star (Burgay 
et al. 2006), that are usually associated with these systems 



(e.g. PSR B1957-F20, Fruchter et al. 1988). Pulsed gamma- 
ray emission seems to come mostly from the 1-3 GeV band 
(Fig. 1). 

In X-rays, the position of PSR J0610-2100 has been 
imaged with the Swift XRT instrument in six observations 
performed in March 2010, April and May 2011. The eflcective 
exposures of those observations were very short, between 0.4 
and 4.1 ks, and with the pulsar position offset from the XRT 
aim point by angles between 0'.2 and 4.'6. PSR J0610-2100 
was not detected with Swift. We used the level 2 data taken 
in Photon Counting mode from these six data sets (10.0 ks 
of the total effective exposure) and applied the approach de- 
veloped by Weisskopf ct al. (2007) for statistical estimates 
on source detections to put a 3 a limit on the XRT source 
count rate of 1.3cntks~^ (in a 30"-radius aperture centred 
at the pulsar radio position). Assuming, for simplicity, a 
power-law X-ray spectrum with a photon index Fx = 2 and 
absorbing hydrogen column density A^h = 0.9 x lO'^^ cm~^ 
(equal to the total Galactic Hi column in the pulsar direc- 
tion^), this count rate translates into an 3cr upper limit of 

1.1 X 10^^ (d/4.0kpc)^ergs"^ for the X-ray luminosity of 
the pulsar in the band 0.3-10 keV. This limit would be more 
stringent if the distance is over-estimated (see the end of 
section 2.2.1). 

PSR J0610-2100 has a 1.4-GHz pulse profile which con- 
sists of 3 components: a main pulse and precursor separated 
by 0.1 rotations and another component which lags the main 
pulse by about 0.3 rotations (see Fig. 1). 

3.2 PSR J1024-0719 

This is a 5.2 ms isolated pulsar located less than 500 pc from 
the Sun (Bailes et al. 1997). An X-ray counterpart was pro- 
posed by Becker & Triimpcr (1999) based on _R05ylr obser- 
vations and later, using XMM-Newton data, Zavlin (2006) 
detected pulsed emission. A candidate optical counterpart 
was reported by Sutaria et al. (2003) but the association 
could be due to positional coincidence and further observa- 
tions arc necessary. 

The gamma-ray pulse profile of PSR J1024-0719 ex- 
hibits one broad gamma^ray peak, roughly half a rotation 
wide, with its centre preceding the radio pealc by about 0.5 
rotation. Pulsed emission seems to be present mostly in the 
1-3 GeV band. 

We include the X-ray pulse profile of PSR J1024-0719 
in Fig. 2 obtained from data collected with the EPIC-pn 
instrument, operated in Timing mode, in a XMM-Newton 
observation conducted in December 2003 for a 66-ks effec- 
tive exposure (Zavlin 2006). We processed the data as in 
Zavlin (2006) but using the latest XMM-Newton data re- 
duction software (SAS v. 11.0.0). The light curve was ob- 
tained from 815 events extracted from columns 38-39 in the 
one-dimensional EPIC-pn CCD image and in the 0.3-2 keV 
energy range to maximise the signal-to-noise ratio. We used 
the photons plugin^ for tempo2 to assign a phase to each 
selected photon, together with a timing ephemeris obtained 
from JBO observations between November 2002 and Octo- 
ber 2004. 

^ http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
^ http://www.physics.mcgill.ca/~aaj:chiba/photons_plug.html 
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Table 3. Gamma-ray light curve and spectral properties for the 6 detected MSPs. 
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Parameter J0610-2100<^ J1024-0719 J1600-3053 



Weighted i?-test statistics 


50.872 


45.405 


103.794 


Pulsation significance ((t) 


6.05 


5.69 


8.84 


Gamma-ray peak multiplicity 


1 


1 


1 


Radio-to-gamma-ray lag, 5 


0.57 ± 0.01 


0.6 ± 0.2 


0.17 ± 0.02 


Gamma-ray peak separation, A 








Photon index, V 


1.2 ± 0.4 

— 0.1 


1.3t 


1.3t 


Cutoff energy, Ec (GeV) 


1.6 ± 0.8 to 2 


2.2 ± 0.7 +'^1 


5.0 ± 1.8 tfi" 


Photon flux, F 0.1 GeV, 10-^ ph cm-2 s"!) 


7.8 ± 2.5 +!!'^ 


4.0 ± 1.1 til 


3.5 ± 0.9 tf°3 


Energy flux, G (5= 0.1 GeV, IQ-^^ erg cm-2 s"!) 


6.6 ± 1.1 


3.8 ± 0.7 


5.2 ± 1.0 tol 


Luminosity, L^,//s2 (J: 0.1 GeV, 10^3 erg s"!) 


10 ± 6 te 


0.11 ± 0.05 +?;!',;' 


4 ± 3 


Efficiency, ri/fa 0> 0.1 GeV) 


11 ± 27 til 


>0.8 ± 0.3 to't 


0.5 ± 0.4 i";* 


Parameter 


J1713+0747 


J174H-1351 


,12051-0827 


Weighted fl-test statistics 


53.675 


51.662 


51.718 


Pulsation significance (cr) 


6.23 


6.10 


6.10 


Gamma-ray peak multiplicity 


1 


1 


1 


Radio-to-gamma-ray lag, <5 


0.32 ± 0.05 


0.74 ± 0.01 


0.54 ± 0.04 


Gamma-ray peak separation, A 








Photon index, P 


1.6 ± 0.3 tal 


1.3t 


1.3t 


Cutoff energy, _Ec (GeV) 


2.7 ± 1.2 


3.1 ± 1.6 li^3^ 


2.0 ± 0.5 to,t 
4.3 ± 1.0 l^J 


Photon flux, F (j! 0.1 GcV, lO"'^ ph cm-^ g-l) 


13.3 ± 3.7 2 


2.8 ± 1.1 til 


Energy flux, G 0.1 GeV, IQ-i^ erg cm-^ s"!) 


10.2 ± 1.5 


3.1 ± 0.8 t°oi 


3.8 ± 0.7 to.t 


Luminosity, L^/fn 0.1 GeV, 10^3 erg s'^) 


1.5 ± 0.3 toi 


0.43 ± 0.12 torn 


0.5 ± 0.2 to^ 


Efficiency, r)//n 0.1 GeV) 


0.44 ± 0.10 


0.020 ± 0.005 toxml 


0.10 ± 0.03 to^l 



Note. — The first quoted uncertainties are statistical, while the second are systematic. Values marked with the symbol 
t were fixed in the spectral analysis and in these cases the systematic errors were calculated in a different manner. 
Details on the measurement of these parameters are given in section 2.2. For PSR J1024— 0719, the reported gamma-ray 
efficiency r; is a lower limit, since only an upper limit on the Shklovskii-corrected period derivative, P, is known (sections 
2.1.1 and 2.2.1). ''Please see the discussion about this pulsar's efficiency in section 2.2.1. 



The H-test reports a, 7.5 a detection of the pulsed emis- 
sion, confirming the result of Zavlin (2006). The estimated 
intrinsic pulsed fraction, corrected for the background con- 
tribution, is 54 ± 20%. The X-ray pulse peaks at phase 0.1- 
0.15, as determined by fitting a series of harmonics to the 
pulse profile, and its shape indicates a predominantly ther- 
mal origin of the pulsar X-ray emission (see Zavlin 2006, for 
more details). 

At ~ 1.4 GHz, this pulsar presents a complex average 
radio pulse profile consisting of about 9 components span- 
ning 0.2 rotations, which emerge from a wider base almost 
half of a rotation wide. The profile at 2 GHz looks very sim- 
ilar (Fig. 2). The emission at 1.4 GHz is almost completely 
linearly polarised in the leading part of the pulse, involving 
the three main peaks, but it is not polarised on the trailing 
part (Yan ct al. 2011). There is little position angle (PA) 
variation detected across the profile. 



3.3 PSR J1600 3053 

This is a 3.6 ms pulsar in a 14.4 day orbit with its binary 
companion (Jacoby et al. 2007). Its gamma-ray peak is ~ 0.3 
rotations wide and exhibits a sharp leading edge and a slowly 
decaying trailing edge. Most pulsed emission comes from the 
two higher energy bands, although some pulsed emission at 
earlier phases may also come from the 0.1-1 GeV band. 

PSR J1600-3053 was observed with XMM-ATewton in 
February 2008 with the EPIC-MOS and EPIC-pn instru- 



ments, operated in Full Window and Timing mode (respec- 
tively), for ~ 30-ks effective exposures. Data were reduced 
with the SAS, version 11.0.0. Examining the EPIC-MOS im- 
ages around the radio position yielded no detection of the 
pulsar. The emission detected in the one-dimensional EPIC- 
pn image was heavily contaminated by enhanced back- 
ground and by another bright source in the field of view. 
Therefore, despite the much higher sensitivity of the EPIC- 
pn instrument, no timing analysis of these data turned out 
to be meaningful (see the example of PSR J0034— 0534 in 
Zavlin 2006). Using the approach described in section 3.1, 
we put a 3 a upper limit of 1.2cntks~^ on the pulsar EPIC- 
MOS count rates, in the 0.3-10 keV band (as measured from 
a 30"-radius aperture centred at the pulsar radio position). 
For the power-law X-ray model of Fx = 2 and absorb- 
ing hydrogen column density A''h = 1.0 X 10^^ cm^^, these 
count rates translate into a rather deep 3 a upper limit of 
7.1 X 10^° (d/2kpc)^ ergs"^ for the X-ray luminosity, in the 
0.3-10 keV band. 



The average pulse profile at 1.4 GHz is dominated by 
a sharp pulse which is preceded by a broader additional 
component, less than 0.1 rotations earlier, with roughly half 
its amplitude. There are no major differences with the pulse 
profile at 2 GHz (Fig. 3). The emission at 1.4 GHz is 30% 
linearly polarised and exhibits two orthogonal PA jumps, the 
second one coincident with a sense reversal of the circular 
polarisation (Yan et al. 2011). 
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Figure 1. Gamma-ray and radio pulse profiles for PSR 
J0610— 2100. The profiles are shown twice, here and elsewhere 
in the paper, to elucidate their complex nature and to clearly 
visualise alignment or otherwise. The gamma-ray light curve for 
events with energies above 0.1 GeV is shown in the top panel. 
The following panels show the light curve for different gamma- 
ray energy bands, as indicated. The histograms were made using 
photon weigfits, and the horizontal dashed lines represent the es- 
timated background emission (section 2). The bottom plot shows 
the radio pulse profile at the indicated frequencies. 
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Figure 2. Gamma-ray and radio pulse profile for PSR 
J1024— 0719 (see Fig. 1 for details). The X-ray light curve is pre- 
sented on the fifth panel from top to bottom. 



3.4 PSR J1713-I-0747 

This is a 4.6 ms pulsar in a 68 day orbit with a white dwarf 
(Foster et al. 1993). The gamma-ray pulse profile of PSR 
J1713-I-0747 is wide and roughly triangular but not exactly 
symmetric, with a shallower trailing edge (Fig. 4). 

In the X-ray band, the position of PSR J1713+0747 has 
only been observed in four short exposures by Swift/XBT 
in Photon Counting mode (3.1 ks of the total effective 
exposure). Using these XRT level 2 data and the same 
approach as in section 3.1 yielded a 3cr limit of 1.9 x 
lO'"*^ (rf/l.lkpc)2 ergs ^ on the X-ray luminosity of this 
MSP in 0.3-10 keV (adopting A^h = 0.5 x 10^^ cm'^). 

The average pulse profile of PSR J1713-I-0747 at 

1.5 GHz consists of one sharp pulse with a weak trailing 



component and two small components on its base, preced- 
ing the main pulse by 0.04 and 0.08 rotations (Yan et al. 
2011). Beside some broadening at lower frequencies, there 
are no major differences between the pulse profile at 0.8 and 
2 GHz. There is a third, very shallow component visible in 
the three radio bands trailing the main peak by a little more 
than 0.1 rotations (sec Fig. 4). The emission at 1.4 GHz is 
almost 100% linearly polarised at the leading and trailing 
edges of the profile. There are two orthogonal PA jumps at 
each side of the main peak and a third one just before the 
shallow trailing component. The second jump (and possi- 
bly also the first) is coincident with a sense reversal of the 
circular polarisation (Yan et al. 2011). 
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Figure 3. Gamma-ray and radio pulse profile for PSR 
J1600-3053 (see Fig. 1 for details). 
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Figure 4. Gamma-ray and radio pulse profile for PSR 
J1713+0747 (see Fig. 1 for details). 



3.5 PSR J1741-I-1351 

This is a 3.8 ms pulsar in a 16 day period binary sys- 
tem (Jacoby et al. 2007). The gamma-ray peak of PSR 
J1741+1351 spans ~ 0.3 rotation and exhibits relatively 
sharp edges. Most of the emission is in the two upper energy 
bands and leads the radio peak, an unusual situation among 
gamma-ray pulsars. Other cases like this arc the MSPs PSRs 
J1744-1134 and J2124-3358 (e.g. Abdo et al. 2013). 

PSR J1741 + 1351 was in the field of view of the Swift 
XRT instrument in January 2012 for 3.7 ks. However, as it 
was projected very close to the edge of the XRT image (with 
a 7.' 3 offset angle), no meaningful limit for the pulsar X-ray 
flux could be derived. 

The pulse profile of PSR J1741-hl351 at 1.4 GHz ex- 
hibits one sharp main pulse and a smaller 0.1 rotation wide 
additional component, leading the main pulse by ~ 0.35 ro- 
tations (see Fig. 5) . There is also a small component on the 
base of the main pulse, which appears more prominent in the 
0.3 GHz profile. This component precedes the main pulse by 
~ 0.05 rotations. 



3.6 PSR J2051-0827 

This 4.5 ms pulsar is one of the so-called Black Widow sys- 
tems. It is in a 2.4 hour orbit and exhibits radio eclipses for 
about 10% of the orbital period, observable at low frequen- 
cies (< 0.6 GHz, Stappers et al. 1996). 

The gamma-ray profile of PSR J2051-0827 exhibits one 
steep, ^ 0.4 rotation wide, peak preceding the radio peak 
by half a period. The pulsed emission is similar in the three 
energy bands (Fig. 6). This profile is compatible with the 
light curve presented by Wu et al. (2012). 

XMM-Newton observed PSR J2051-0827 in April 2009 
with the EPIC-MOS and EPIC-pn detectors, all operated in 
Full Window mode, for 44 and 37-ks effective exposures, re- 
spectively. Unfortunately, these observations suffered from 
numerous strong particle fiares that increased the back- 
ground level by a factor of up to 10-12 (compared to the nor- 
mal one) during about 80% of the observations. This made 
the XMM- iVeuiton data practically useless for a purposeful 
analysis. 

Five observations of this object were conducted with the 
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pulsar dispersion measure. Fixing Nh ~ 0.6 x lO'^^ cm ^ 
results in a power-law model with Fx — 2.6 and unab- 
sorbed flux Fx ~ 7.9 x 10"^^ ergcm"^ s"^ in the 0.3-10 keV 
range, in agreement with the result reported by Wu et al. 
(2012). The thermal blackbody model provided more rea- 
sonable parameters; Nn < 1 x lO^'^cm"^, a (redshifted) 
temperature ~ 2.7 x 10® K and radius of the emit- 
ting area (hot polar caps on the pulsar's surface) ~ 
0.04 (d/l kpc) km. The corresponding bolometric X-ray lu- 
minosity is LjJ ~ 0.6 x 10"'"' (d/l kpc)^ crgs^"*^. Using a non- 
magnetic hydrogen atmosphere model for thermal emission 
of MSPs (Zavlin et al. 1996; Zavlin 2006, 2009), we ob- 
tained Tpc ~ 1.6 X 10*^ K, Rpc ~ 0.18(fi/lkpc)km and 
Lpc ~ 0.8 X 10^" (d/lkpcf crgs"\ Tlicsc arc all uincd- 
shifted values, i.e. as measured at the neutron star surface. 

Using the photons plugin for tempo2, we searched for 
orbital variability of the X-ray photons detected in the 
Chandra observations. With the scant count statistics avail- 
able, we found that the significance of a possible modulation 
does not exceed 2.2 a. X-ray data of much better quality are 
required to elucidate the origin of the emission (whether it 
is thermal polar cap radiation or non-thermal flux from the 
pulsar magnetosphcrc and/or interaction of the pulsar wind 
with the companion). 

The pulse profile of this MSP at ~ 1.4 GHz exhibits 
one main pulse and one overlapping trailing component 
about 0.05 rotations apart. Additionally, there are at least 
3 smaller trailing components extending up to 0.2 rotations 
away, forming a shoulder. The pulse profile at 0.4 GHz keeps 
the same principal structure but the two main peaks appear 
broader and slightly more separated. There are no major dif- 
ferences with the profile at 2 GHz (sec Fig. 6). The emission 
at 1.4 GHz is mildly linearly polarised (~12%) throughout 
the whole profile (Xilouris et al. 1998). 



Figure 5. Gamma-ray and radio pulse profile for PSR 
J1741-t-1351 (see Fig. 1 for details). 



Chandra ACIS-S instrument in Very Faint mode in March 
and July 2009 for 9-ks effective exposures (see also Wu et al. 
2012). The aim of these observations was to search for or- 
bital variability of the pulsar X-ray flux and did not pro- 
vide a time resolution suitable to perform a timing anal- 
ysis at the pulsar spin period. We reduced the data using 
the CIAO software, version 4.3. Despite the small number 
of source counts collected in each observation (8-12 counts 
in a l"-radius aperture centred at the pulsar radio posi- 
tion), the object was clearly detected. The combined dataset 
totalled 44 source counts in the 0.3-10 keV range, with a 
negligible (less than 1%) background contamination. The 
spatial distribution of these counts was found to be consis- 
tent with the ACIS-S point-like source image. Of the ex- 
tracted source counts, 90% were detected at photon en- 
ergies below 2 keV, indicating that the spectrum of PSR 
J2051— 0827 is soft and likely of a thermal origin rather 
than of a non-thermal (magnctospheric) one. Indeed, fit- 
ting a power-law model to the extracted spectrum resulted 
in a large best-fitting photon index Fx — ^ and a hy- 
drogen column density JVh ^ 2.3 x 10^^ cm~^ significantly 
greater than the value of 0.6 x 10^^ cm~^ expected from the 



4 SEARCH FOR ORBITAL MODULATION OF 
THE GAMMA-RAY EMISSION 

PSRs J0610-2100 and J2051-0827 are in close binary or- 
bits, where interactions between the wind of the pulsar and 
the atmosphere of the companion star are likely to occur. 
Accordingly, we have searched for gammas-ray flux modula- 
tions as these systems follow their orbital cycles. For com- 
pleteness, we also searched for such modulations in the other 
three binary systems, PSRs J1600-3053, J1713-I-0747 and 
J1741-I-1351. 

The Fermi-LAT sensitivity to a source varies on multi- 
ple timcscales, most notably the spacecraft's orbital period 
(~95 minutes) and period of precession (~53 days). Beats 
between the binary orbital frequency and harmonics of these 
timescales induce an apparent modulation of photon rate as 
a function of orbital phase. Searches for intrinsic modulation 
must therefore carefully correct for the time-dependent sen- 
sitivity (e.g. Corbet & Kerr 2010; Guillemot et al. 2012b). 

Because of the relatively low count rates, wc decided to 
adapt an unbinned pulsation statistic, the H-test (de Jager 
et al. 1989), for use with uneven exposure. To do this, we 
computed the exposure to the source with 30 s resolution 
over broad energy bands (2 per decade). We used Tempo2 
and the timing solutions of §2.1.1 to compute the orbital 
phase for each 30 s interval and took the resulting distri- 



New Fermi MSPs and the radio /gamma-ray connection 11 



16 

14 

12 

<5 10 
c 

g 8 
o 

^ 6 
4 

2 

2.5 
1.5 




o 
o 



1 T 

0.5 y 



i2 6 
g 4 



S 2 -J 



§ 6 
8 4 



S 2 



H — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 



> 3 GeV 




0.1-1 GeV 



L NRT@2GHz JV / 




GC 



0.5 1 1.5 2 

Pulse phase 

Figure 6. Gamma-ray and radio pulse profile for PSR 
J2051-0827 (see Fig. 1 for details). 



bution of phases, J^, to represent the null hypothesis of no 
intrinsic modulation. By definition, the quantity with 
(j) the observed orbital phases, is a uniformly distributed 
random variable in the absence of intrinsic modulation, and 
thus suitable for use in the H-test. 

Using this modified H-test, along with the photon 
weights (see §2.2), we searched for a signal in individual 
energy bands (100-300, 300-1000, 1000-3000, 3000-10000, 
and 10000-30000 MeV) as well as cumulatively (>100, >300, 
<300 MeV, etc.). We detected no significant orbital modu- 
lation of the gamma-ray signal from any of the sources. 

This null result is physically expected for the binary sys- 
tems with long orbits, where no direct interaction between 
the two bodies is feasible. In the case of PSR J2051-0827, 
this result is also unsurprising as no evidence for interac- 
tion of the pulsar wind with the companion star is seen 
in the X-ray spectrum (sec section 3.6). Interestingly, the 
weighted orbital light curve at energies above 3 GeV for PSR 
J0610— 2100 presents two peaks, suggesting the presence of 
orbital modulation of the emission at these energies (Fig. 




Orbital phase 

Figure 7. Gamma-ray emission above 3 GeV folded on the orbital 
period for PSR J0610— 2100. The light curve was made using 
photon weights and it is shown twice for clarity. Vertical dotted 
lines indicate the phase at which eclipses would be expected and 
the horizontal dashed line represents the photon background. The 
H-test for this light curve is low, suggesting the modulation is not 
significant. 



7). However, there are very few events above 3 GeV and the 

peaks observed could be the effect of low statistics; indeed, 
the H-test for this set of events indicates only a 2 a signifi- 
cance. Furthermore, the light curve of the background pho- 
tons (obtained by assigning a weight 1 — to every event in 
the same 5 deg ROI, see section 2.2.2.) for this source seems 
to follow the same shape in Fig. 7. Although fairly flat, the 
exposure folded curve for PSR J0610— 2100 also presents lo- 
cal peaks at similar orbital phases (~ 0.2 and ~ 0.5). 

Finally, we performed a Monte Carlo study to deter- 
mine the sensitivity to a few particular types of modulation. 
For each pulsar, we use the observed weights and an as- 
sumed morphology for the orbital modulation to generate 
100 random realisations of the orbital phases. We then vary 
the strength of the modulation until 95 of the 100 simula- 
tions exceed the statistical threshold, taken to be H = 8, or 
a ~2(T detection. 

We find in general that only strong modulations are 
detectable. Sinusoidal modulation is detectable for PSRs 
J0610-2100 and J1713-h0747 if ~70% of the total fiux is 
modulated, while not even 100% sinusoidal modulation is 
reliably detectable for the remaining sources. On the other 
hand, such high levels of modulation, which imply addi- 
tional, non-magnetospheric sources of gamma-ray emission 
from the system, are incompatible with the estimated back- 
ground levels of the rotational phase light curves (Figs. 1 to 
6). 

For the case of a notch in otherwise-steady emission 
(an eclipse), we find that the sensitivity depends strongly 
on the notch width. Any gamma-ray eclipse is expected to 
be narrow, and unfortunately the LAT is largely insensitive 
to such small features. Eclipses become detectable for PSRs 
J0610-2100 and J1713-F0747 when they span ~20% of the 
orbit, while even greater values (30-50%) are required for the 
remaining sources. Thus, while we can rule out extreme cases 
of modulation (fully-modulated sinusoids, broad eclipses), 
our sensitivity study shows the null results discussed above 
are not highly constraining. 
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5 RADIO AND GAMMA-RAY PROPERTIES 
OF THE GAMMA-RAY MSP POPULATION 

Our sample of 30 MSPs is constituted by all MSPs known 
to emit in gamma rays at the time of writing (published 
or publicly announced). We have compiled published radio 
pulse profiles and measurements of the radio spectral in- 
dex (s), radio flux densities, average degree of radio linear 
polarisation through the main radio pulse, phase lag of the 
closest gamma-ray peak to the main radio peak (6), gamma- 
ray photon index (F) and flux for each of the 30 MSPs in the 
sample. Table 4 gives some of these quantities and also lists 
references where gamma-ray and radio pulse profile plots 
can be found. 



5.1 Radio and gamma-ray pulse profile properties 

Considering properties like phase-alignment between radio 
and garnrna-ray peaks (above 0.1 GcV), radio duty cycle 
and number of pulse components, we define three main 
groups and label them A, N and W (Table 4). Although 
this analysis is based mainly on the inspection of average 
pulse profiles at ~ 1.4 GHz, we also inspected pulse profiles 
at lower and higher radio frequencies. None of the MSPs 
studied exhibits appearance or disappearance of pulse com- 
ponents (though, see comments on PSR B1957-I-20 below 
and on PSR J0218+4232 in section 5.1.2). We note that 
pulse-profile evolution with frequency is observed in some 
other MSPs. 



5.1.1 Type A(ligned) 

These are the MSPs which have their main gamma-ray peak 
aligned with the main radio pulse. There are five objects in 
this group (Table 4). 

All of the A-type MSPs exhibit two main radio peaks, 
with each of them composed of one or more components. 
In two cases (PSRs B1937-t-21 and B1957-f20) the second 
peak appears as an interpulsc, i.e. half a rotation away from 
the main pulse. However, in the case of PSR B1937-I-21, low 
level emission detected between the peaks suggested that 
both peaks may not be due to emission produced at op- 
posite magnetic poles but somewhere in the outer magne- 
tospherc (Yan et al. 2011). Only five MSPs in the whole 
sample present secondary peaks about half a rotation away 
from the main pulse. 

For the five A-type MSPs, all radio peaks have a 
gamma-ray counterpart. The only exception might be PSR 
B1957-I-20, for which the second radio peak visible at 
1.4 GHz, leading the main peak by ~ 0.26 rotation, appears 
not to have an obvious gamma-ray counterpart. We note 
that this peak is not visible at lower radio frequencies (e.g. 
0.3 GHz). 

Four MSPs in this group exhibit very low or unde- 
tectable levels of linear polarisation in their radio emission 
(Table 4). This seems to be a property exclusive to A-type 
MSPs: all other gamma-ray MSPs in the sample for which 
polarisation data were available exhibit normal levels of lin- 
ear polarisation. The only A-type MSP showing some degree 
of linear polarisation is PSR B1937-I-21. 



5.1.2 Type N(on-ahgned) 

These are MSPs having their main gamma-ray peak out 
of phase with the main radio peak and their radio emission 
dominated by a single peak. Two MSPs in this group present 
gamma-ray pulses almost phase-coincident with their main 
radio pulses (PSRs J1744-1134 and J2214+3000). However, 
they are not exactly aligned and their peaks lead the radio 
peak, which is unusual among gamma-ray pulsars; though it 
is also seen in the fight curves of the MSPs PSR J1741-M351 
and PSR J2124— 3358. These objects may belong to a differ- 
ent class of MSPs, in terms of morphology (cf. Venter et al. 
2009; Johnson 2011). 

There seems to be some gamma-ray emission, besides 
the main peak of emission, at the radio pulse phase in 
at least 3 objects (PSRs J0340-I-4130, J0610-2100 and 
,11024-0719). There are 21 N-typo MSPs in total, two 
of them exhibiting an interpulse (PSRs J0030-I-0451 and 
JOlOl-6422). 

5.1.3 Type W(ide) 

This is a group of four MSPs that have their main gamma- 
ray peak out of phase with the main radio peak and 
whoso radio emission consist of multiple peaks of compa- 
rable amplitude, covering most of the rotational period. 
PSR J0218-I-4232 presents a secondary radio pulse com- 
ponent coincident with the gamma-ray peak of emission. 
This radio component is stronger at lower radio frequen- 
cies (< 0.6 GHz), becoming comparable in strength to the 
main radio peak (Kuiper et al. 2002). There also seems to 
be alignment between some secondary radio components 
and the main gamma-ray peak for PSRs J2124— 3358 and 
J2302-I-4442. Nevertheless, they are not proper A-type pul- 
sars because not all of their radio peaks are aligned with 
gamma-ray peaks. In addition, it must be noted that given 
the large number of radio components and their long duty 
cycles, those could be mere coincidences. Besides these prac- 
tical reasons, there is a physical reason to keep the W and N 
groujjs separated. The W-typo pulsars are possibly aligned 
rotators, i.e. pulsars with their magnetic axes almost coin- 
cident with their rotation axes, a situation that could affect 
their observed properties. 

5.2 Magnetic field at the light cylinder 

The dipolc magnetic field strength of a pulsar is normally 
calculated by assuming that all rotational energy losses are 
due to dipole radiation. Its value at the light cylinder can 
be estimated by (e.g. Lorimer & Kramer 2005) 

At the light cylinder, the co-rotational velocity is equal 
to the speed of light, c, and its radius is therefore given by 
Rlc = cP/27r, where P is the rotation period. This implies 
that gamma-ray MSPs have the smallest light cylinder radii, 
followed by rapidly spinning young pulsars like the Crab 
pulsar. Because of the weak dependance of Blc on P and 
the different P ranges exhibited by MSPs and young pul- 
sars, both populations share similar Blc values, the largest 
among the whole pulsar population (> 10* G). 
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Table 4. Radio and gamma-ray properties of 30 MSPs detected by the LAT. 
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Note. — Columns contain: pulsar name, magnetic field at the light cylinder (S^c), mean degree 
of linear polarisation {L/I), radio spectral index s, radio lag (<5) and type of MSP, according 
to the classification in section 5. Uncertainties on the last quoted digit are indicated between 
parentheses. For the average linear polarisation, the errors quoted correspond to a rough and 
conservative estimate of the standard deviation from several measurements made by different 
authors, preferably around 1.4 GHz. Type labels marked with the symbol f are for those MSPs with 
radio and gamma-ray pulses overlapping in phase but not considered as aligned (section 5.1.1). The 
4 references in the last column are for {L/I), s, 5 and for a gamma-ray/radio pulse profile plot: (1) 
Lommen et al. (2000); (2) Stairs et al. (1999); (3) Kerr et al. (2012); (4) Navarro et al. (1997); (5) 
Yan et al. (2011); (6) Ord et al. (2004); (7) Manchester & Han (2004); (8) Xilouris et al. (1998); 
(9) Keith et al. (2012); (10) Camilo et al. (2012); (11) Thorsett & Stinebring (1990); (12) Pruchter 
et al. (1990); (13) P. Demorest, S. Ramson (priv. communication) (14) Toscano et al. (1998); (15) 
Navarro et al. (1995); (16) Bangalc et al. (in preparation) (17) S. Ransom (priv. communication) 
(18) Kramer et al. (1998); (19) Demorest et al. (2012); (20) M. Keith (priv. communication) (21) 
Abdo et al. (2013); (22) Abdo et al. (2010a); (23) Keith et al. (2011); (24) Guillemot et al. (2012b); 
(25) Abdo et al. (2010b) (26) Ransom et al. (2011) (27) Preire et al. (2011a) (28) Cognard et al. 
(2011) (29) Guillemot et al. (2012a) (30) This work; 



A-type MSPs tend to have higher Blc values than the 
rest of the population. Fig. 8 shows B^c as a function of 
P for all known pulsars having P ^ 0.1s. If proper mo- 
tion measurements were available through the ATNF Pulsar 
Catalogue, the Blc values were corrected for the Shklovskii 
effect (Table 4). Lines of constant spin-down energy rate 
E = Att^IP/P^ (where / is the moment of inertia of the star, 
assumed to be 10'*®gcm^) are plotted with dashed lines. 
Most MSPs populate the lower branch crossing the centre 
of the plot and the 30 gamma-ray MSPs in our sample pop- 
ulate its higher end, towards the shortest periods, and larger 
E and Blc values. The MSPs from the three types described 



in the last section are plotted using different symbols. A-type 
MSPs are near the top left corner, with the shorter periods 
and the higher Blc and E values (Johnson 2011). 

PSR B1823-30A is one of the A-type MSPs in the plot 
and, considering its period, shows a relatively large Blc- 
This MSP is in a globular cluster (GC) and has a very large 
P, hence its large Blc- Freire et al. (2011a) studied the 
gamma-ray emission of this MSP and concluded that the 
observed P is intrinsic to the pulsar and not due to acceler- 
ation in the cluster. MSPs in GCs are marked in Fig. 8 and, 
due to possible contamination of their P values produced by 
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Figure 8. Magnetic field at the light cylinder, corrected for the Shklovskii effect if possible, versus the rotational period for all known 
pulsars with P 0.1s. MSPs in globular clusters (GC) are shown using thick black dots to indicate that their Blc values are uncertain 
(see section 5.2). Lines of constant spin-down energy rate are dashed. 



movement in the gravitational potential of a GC, their Blc 
values are uncertain. 

A 2-dimensional Kolmogorov-Smirnov (KS) test 
(Fasano & Franceschini 1987; Press et al. 1992) be- 
tween A-type MSPs and N and W types combined 
indicates that the probability that they belong to the 
same P-Blc distribution is less than 1%. However, the 
significance of the 2-dimcnsional KS test remains valid 
while = iViiV2/(A''i + iV2) ^ 20 (Press et al. 1992) and, 
in this case, A^e = 4.2. To check these results further, we 
also tried the standard 1-dimensional KS test, that can 
be applied for lower number of objects (valid for ^ 4, 
Press et al. 1992). Comparing the BLC-values of the five 
A-type MSPs with the 25 N and W types gives 0.02% 
probability that they belong to the same distribution. To 
quantify how likely it is to obtain such a value by selecting 
a small group at random, we randomly picked five MSPs 
from the sample and calculated the KS test against the 
remaining 25 objects. After repeating this 30,000 times 
we obtained that for more than 99% of the cases the KS 
test null hypothesis probability is greater than 0.02%, with 
67% of the cases having a probability greater than 40%. 
Therefore, we conclude that the distribution of BLC-values 
of the A-type MSPs is significantly difl^erent to that of the 
rest of the gamma-ray MSPs. 

We acknowledge, nonetheless, that the above results 
might be driven by the very short periods that most of 
the A-type MSPs have, compared to those of the N and W 



types together. Indeed, while the KS test gives a probability 
of 41% for the P values being from the same distribution, 
the same test gives only 0.4% probability that the periods 
belong to the same distribution. However, this is about 10 
times larger than the probability obtained for Blc- If E val- 
ues are considered instead, the probability is 1%, which is 
again 50 times larger than what is obtained for Blc values. 
Moreover, as will be discussed in section 6, the comparison 
with the Crab pulsar plus other considerations hint towards 
Blc being the relevant parameter, rather than P or E. 

5.3 Radio spectral indices 

The radio spectral index s corresponds to the slope of the 
spectrum in a logarithmic scale, in which the flux density is 
described hy S oc j/" . Out of the 30 MSPs in the sample, we 
have spectral information for 19 of them (Table 4). These 
data show that A-type MSPs tend to have lower spectral 
indices than the rest of our sample. Below we assess the 
significance of this observation. 

The distribution of spectral indices for all the gamma- 
ray MSPs together with 33 additional MSPs for which spec- 
tral information was available is shown in Fig. 9, indicat- 
ing the contribution made by the A-type MSPs. Data were 
collected from the ATNF Pulsar Catalogue, from the com- 
pilations in Lorimer et al. (1995), Kramer et al. (1998) and 
Toscano et al. (1998), and from flux densities at different fre- 
quencies obtained via private communication with M. Keith 
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Figure 9. Distribution of radio spectral indices for 52 MSPs 
highligliting the different values for gamma-ray MSPs and A-type 
MSPs. 



and S. Ransom. Note that the error bars of the s-values 
sometimes may not reflect the intrinsic difficulties associated 
with flux density measurements, which directly affect the 
spectral index estimate. While the s distribution of gamma- 
ray MSPs seems to follow the general distribution, with a 
mean spectral index above —2.0 (consistent with Kramer 
et al. 1998; Toscano et al. 1998), the mean value for the 
A-type gamma-ray MSPs clearly falls below —2.5 (see Fig. 
9). 

To test how significant the difference between the s dis- 
tributions of the different sets is, we again use the KS tost. 
A KS test over the s values of four A-type MSPs against 
the remaining 15 gammarray MSPs (Table 4) indicates a 
probability of 0.6% that both sets belong to the same distri- 
bution. If the test is applied for the A-type MSPs against all 
other MSPs in Fig. 9 the probability falls to 0.2%. There- 
fore, it appears that the A-type MSPs have distinctively low 
s-values. We note, however, that for the first KS test men- 
tioned above, A^e = 3.2 (see section 5.2). On the other hand, 
the test indicates a probability of 56.7% that the s values 
of all gamma-ray MSPs belong to the overall s distribution, 
indicating that in general there is no evidence for different 
radio spectral behaviour between gamma-ray MSPs and the 
rest of the MSP population. 

We also checked how likely it is to obtain an s distri- 
bution such as the one exhibited by the A-type MSPs by 
pure chance. To do so, we randomly picked four values from 
the s-distribution of all gamma-ray MSPs and calculated 
the probability of them being from the same distribution as 
the remaining 15. After repeating this process many times, 
we find that more than 99% of the trials give a probabil- 
ity greater than 0.6%, with more than 66% of the trials 
giving a probability greater than 40%. Thus, the particular 
s-distribution of A-type MSPs appears to be truly different. 

We find no correlation between s and the gamma-ray 
spectral index or photon index, F, for the MSPs in the sam- 
ple. Nonetheless, we note that values of F are tightly clus- 
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Figure 10. Distribution of the average levels of radio linear po- 
larisation for 18 MSPs in the sample. 

tered, with a scatter smaller than unity (Abdo et al. 2010b). 

As well, no obvious correlation was found between F and any 
other investigated parameter. 

5.4 Fluxes 

There arc radio flux densities at 1.4 GHz available for 23 
of the 30 MSPs in our sample and gamma-ray fluxes for 
all of them. These last values are preliminary results of the 
ongoing effort to produce the second Fermi LAT Catalog of 
gamma-ray pulsars by the Fermi collaboration (Abdo et al. 
2013). For many MSPs in the sample, flux values can also be 
found in the references given for the gammarray light curves 
in Table 4. We find no obvious correlation between radio flux 
densities and gamma-ray fluxes (Ackcrmann et al. 2012b). 

Fluxes are expected to be at legist slightly correlated be- 
cause of their mutual dependance on distance. The fact that 
this correlation is not observed in the sample suggests that 
radio and/or gamma-ray fluxes arc highly dependent upon 
geometrical factors and probably other intrinsic properties 
dependant on P and P. No correlation between these fluxes, 
or their ratio, and any other quantity was found. 

5.5 Linear polarisation 

Wc found published average radio linear polarisation {L/I) 
data for 17 MSPs in our sample. No correlation is found 
between this parameter and any other studied quantity. 
However, we note that the four MSPs with very low (or 
zero) {L/I) are all A-type (PSRs J0034-0534, B1820-30A, 
J1902-5105 and B1957-I-20, see Fig. 10). 



6 DISCUSSION 

6.1 Gamma-ray MSPs at a glance 

By studying the radio and gamma-ray pulse profllcs of 30 
MSPs we have classified gamma-ray MSPs into 3 types, la- 
belled A, N and W. The first group is composed by MSPs 
exhibiting pulsed gamma-ray emission phase-aligned with 
their radio emission. The N and W type MSPs exhibit mis- 
aligned emission, with the N-types presenting radio emission 
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dominated by a single pulse and the W-type presenting mul- 
tiple and wide radio pulses of comparable amplitude. While 
the N and W types appear to constitute the norm among 
gamma-ray MSPs, A-type MSPs appear less frequently and 
present various notable differences from the rest of the pop- 
ulation. We note that these properties are not shared by 
PSRs J1744-1134 and J2214-I-3000, the two MSPs which 
could have possibly been regarded as A-type (see sections 
5.1.1 and 5.1.2). 

W-type MSPs are more likely to be aligned rotators, 
with their rotation axes almost coincident with their mag- 
netic axes. This situation, however, does not appear to drive 
obvious trends either in their gamma-ray profiles or any of 
the other studied parameters. 

The difference in rotational phase between radio and 
gamma-ray pulses has commonly been interpreted as a dif- 
ference in the location, within the pulsar's magnetosphere, 
where the emission is generated. Standard models normally 
place the source of radio emission at lower altitudes com- 
pared to the source of gamma rays, in accordance with the 
radio/gammar-ray misalignment observed for most gamma- 
ray pulsars and MSPs. However, in light of the newly de- 
tected MSPs exhibiting phase-aligned emission, other possi- 
bilities ought to be considered. In the two most likely sce- 
narios, both radio and gammarray emission are generated 
either relatively near the surface of the neutron star (low- 
altitude Slot Gap (laSG) models) or over an extended re- 
gion at higher altitudes (altitude-limited Open Gap (alOG) 
and altitude limited Two Pole Caustic (alTPC) models (cf. 
Venter et al. 2012)). Although the alTPC model seems to 
offer better fits for some A-type MSPs, further studies are 
necessary to better understand the emission mechanisms op- 
erating in these pulsars. 

Nonetheless, models involving a caustic origin for the 
radio emission (at high altitudes) predict low levels of linear 
polarisation and rapid swings of the polarisation position 
angle (Dyks et al. 2004). In such models, rapid position angle 
swings with phase and depolarisation of the total emission 
would be an effect of mixing signals from different altitudes 
(Venter et al. 2011). This is indeed observed in some, but 
not all, of the A-type MSPs (section 5.5; Table 4). 

It has been noted that for high E non-MSP pulsars the 
degree of linear polarisation may be correlated with E (Wel- 
tevrede & Johnston 2008). Moreover, based on their wide 
beams, it was argued that the radio emission of some of these 
pulsars is produced at a large range of heights. Their polar- 
isation levels, however, appear high, contrary to what is ar- 
gued above to explain the lack of polarised emission in some 
A-type MSPs. Also, A-type MSPs have the largest E values 
among MSPs, defying the correlation found by Weltevrede 
& Johnston (2008). Nevertheless, normal pulsars have larger 
light cylinder radii than MSPs and it might be that depo- 
larisation is less efficient in their larger magnetospheres. 

In the smaller magnetospheres of MSPs the production 
of radio emission could naturally extend up to larger frac- 
tions of Rlc, producing the observed wide beams (Yan et al. 
2011) and the phase-aligned profiles of A-type MSPs (Hard- 
ing 2005). In fact, A-type MSPs have shorter rotation pe- 
riods than the other MSPs (Fig. 8). Based on a study of 
the beaming fractions of the radio and gamma-ray emis- 
sion, Ravi et al. (2010) concluded that the radio emission of 
high-£ pulsars (including MSPs) must originate higher up 



in the magnetosphere, near where the gamma-ray emission 
is produced. This is consistent with the high-iJ values of A- 
type MSPs. The Crab pulsar is the only non-MSP known 
to exhibit nearly phase-aligned radio and gamma-ray emis- 
sion (Kuiper et al. 2003). We note that E for the Crab pul- 
sar is almost 3 orders of magnitude larger than the highest 
values among MSPs and that many young gamma-ray pul- 
sars have E values similar to those of the A-type MSPs, 
but show misaligned emission. Also, while the Crab pulsar 
has the secoiid smallest Rlc value aiiioiig normal pulsars, 
it is more than one order of magnitude larger than those of 
the A-type MSPs. If Rlc was the main factor determining 
the alignment between radio and gamma-ray emission, we 
would expect to see this alignment for most MSPs, which 
is not observed. On the other hand, if E was the main fac- 
tor, we would expect to see more young gamma-ray pulsars 
exhibiting aligned radio/gamma-ray emission, which is not 
the case. Considering the KS-tests described in section 5.2 
and the position of the Crab pulsar in Fig. 8, Blc appears 
to be a natural common property among A-type MSPs and 
the Crab pulsar. We note that the polarisation levels of the 
Crab pulsar's radio emission appear to be well above zero 
(Gould & Lyne 1998). 

The Crab pulsar and A-type MSPs all have similar radio 
spectra. The radio spectral index of the Crab pulsar is —3.1 
(Lorimer et al. 1995) , substantially steeper than the average 
of —1.9 for MSPs and of —1.8 for normal pulsars (Toscano 
et al. 1998: Maron et al. 2000). Although viewing angles 
and other geometric factors can bias our measurements, the 
steep radio spectra of the A-type MSPs and the Crab pulsar, 
together with their high Blc values, are likely to be related 
to magnetospheric similitudes and common processes on the 
generation of their emission. 

The emission of Giant radio Pulses (GPs) is another 
common feature among A-type MSPs and the Crab pulsar. 
GPs are sporadic, short and intense bursts of radio emission, 
following power-law energy statistics. They were originally 
detected in the emission of the Crab pulsar and later in 
PSR B1937-F21 (Heiles et al. 1970; Cognard et al. 1996). 
Today there are 8 pulsars and 5 MSPs known to emit these 
type of pulses (Knight et al. 2005; Knight 2006). The CP 
properties exhibited by the Crab pulsar differ from most 
other pulsars but are very similar to those exhibited by the 
5 MSPs found to emit GPs (Knight 2006). Three of these 
MSPs are A-type: PSRs B1937-I-21, B1957-f 20 (Knight et al. 
2006) and B1820-30A (Knight et al. 2005) and the other 
two are PSRs J0218-I-4232 (Joshi et al. 2004), a W-type 
MSP, and B1821-24 (Romani & Johnston 2001), for which 
no gamma-ray pulses have been detected with a confidence 
above 5 a (Pellizzoni et al. 2009). It has been proposed that 
the relatively large Blc values exhibited by these 5 MSPs 
and the Crab pulsar could be the main physical factor de- 
termining the generation and the main properties of their 
GPs (Cognard et al. 1996, but see discussion in Knight et al. 
(2006)). Future observations and careful analyses should dis- 
cern whether the emission of GPs is somewhat connected to 
the alignment between radio and gammarray emission or is 
mere coincidence. No GPs were detected for the A-type PSR 
J0034-0534 (Knight et al. 2005). 

The wide pulse profiles that MSPs generally have, com- 
pared to normal pulsars (Yan et al. 2011), could be un- 
derstood as evidence for outer-magnetosphere caustic radio 
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emission for all types of MSPs. One could argue that the 
outer magnetosphere offers enough room for different emis- 
sion locations for radio and gamma rays. However, we have 
shown that A-type MSPs exhibit different emission proper- 
ties. Could the availability of larger magnetic fields at the 
light cylinder generate conditions, in relatively smaller mag- 
netospheres, that favour the generation of radio emission 
at higher altitudes (at least as a fraction of Rlc) and co- 
located with the production of gamma rays? Would this 
somewhat different emission mechanism naturally produce 
a steeper spectra and favour the production of giant pulses? 
Any model describing phase-aligned radio and gamma^ray 
emission should take these properties into consideration. 



7 SUMMARY 

Wc have presented the detection by the Fermi LAT of 
gamma-ray pulsations from six MSPs, five of them detected 
for the first time and a sixth one, PSR J2051— 0827, con- 
firmed at the 5c7 level. The six MSPs present properties 
which are common among the gamma-ray MSP population. 
All of these pulsars are significantly detected in the radio 
domain, but due to the low conversion of spin-down energy 
into X-rays, only two are detected significantly in X-rays. 

By studying the morphology and phase relationship 
of radio and gamma-ray pulse profiles of a sample of 30 
MSPs, wc grouped gamma-iay MSPs into three types. The 
most distinctive type of gamma-ray MSPs are those exhibit- 
ing phase-aligned radio/gamma-ray emission. We find some 
clear trends in their emission properties, which differ sig- 
nificantly from the rest of the MSP population. We showed 
that the MSPs in this group have a radio spectra steeper 
than the rest of the MSP population and that they also 
have among the highest inferred magnetic field strengths at 
the light cylinder. Additionally, some MSPs in this group 
have distinctively low degrees of radio linear polarisation 
and some of these type are amongst the handful of MSPs 
known to emit giant radio pulses. Many of these properties 
are also observed from the Crab pulsar, the only normal ra- 
dio pulsar known to emit gamma-ray emission phase-aligned 
with its radio emission. 

The use of combined information obtained through the 
study of phase aligned gamma-ray. X-ray and radio emission, 
along with their intrinsic properties, offers a wide perspec- 
tive that certainly helps to improve our understanding of 
the emission mechanism of pulsars. 
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